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ABSTRACT 

Context. This is the third of a series of papers presenting the first attempt to analyze the growth of the bar instability in a consistent 
cosmological scenario. In the previous two articles we explored the role of the cosmology on stellar disks, and the impact of the 
gaseous component on a disk embedded in a cosmological dark matter halo. 

Aims. The aim of this paper is to point out the impact of the star formation on the bar instability inside disks having different gas 
fractions. 

Methods. We perform cosmological simulations of the same disk-to-halo mass systems as in the previous works where the star 
formation was not triggered. We compare the results of the new simulations with the previous ones to investigate the effect of the star 
formation by analysing the morphology of the stellar components, the bar strength, the behaviour of the pattern speed. We follow the 
gas and the central mass concentration during the evolution and their impact on the bar strength. 
Results. In all our cosmological simulations a stellar bar, lasting 10 Gyr, is still living at z=0. 

The central mass concentration of gas and of the new stars has a mild action on the ellipticity of the bar but is not able to destroy it; 
at z=0 the stellar bar strength is enhanced by the star formation. The bar pattern speed is decreasing with the disk evolution. 

Key words, galaxies: spirals, structure, evolution, halos, kinematics and dynamics 



1. Introduction 

The connection between the bar feature and the star for- 
mation process has alre ady been pointed out in the past: 
Martinet & Fr iedlil d 1997b observed that non interacting galaxies 
displaying the highest star forming activity have strong and long 
bars. Conv ersely not all strong and long bars are intensively cre- 
ating stars. iMazzei & CuriH d200ll) performed SPH simulations 
to investigate the dependence of the star formation in disks on 
the geometry and dynamical state of the DM halo. They showed 
that the star formation lengthens the life time of the bar. However 
in their works the evolution of the disk+halo system arises in a 
isolated framework, outside the cosmological scenario. 
The treatment of the star formation in galaxies involves a large 
number of physical processes which arise with different time 
scales and spatial lengths. Therefore i ts inclusions in s imula- 
tions of galaxy formation is complex (M azzei & CuriH (2003) 
and references therein) and many sub-grid parametrisations (i.e. 
representing all the astrophysical processes that cannot be re- 
solved due to resolution or computing power limitations) of the 
star formation process, exist in literature. To mimic all these 
processes with a numerical model of a realistic disk galaxy is 
difficult too. Thus the numerical work on this subject focuses 
either on a detailed analysis of the interstellar medium (ISM) 
and a smaller simulated area dWadal 120011: lAvillez et alJl2000h 
or on the study of the global instabilities and of the star forma- 
tion at the cost of simplified ISM model s (Robertson et al. 2004; 
ISemelin & Combesl 120021: ICombes & Melchiotj|2002l) . Models 
that h ave tackled both topics have either been in two dimen- 
sions dWada & Normanl | 2001l) or re s tricted to a box size a few 
hundred parsecs across dWadall200l|) . |Rosen & Bregmanld 19951) 
worked in two dimensions allowing the ISM to evolve self- 
consistently but neated the stars as collisional rather than col- 
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lisionless fluid. iTasker & Brvanl ((2006) presented the first three 
dimensional simulations of a global disk without the need to sim- 
plify the structure of the ISM. Their paper is devoted to under- 
stand the fundamental processes of star formation and feedback 
in a disk galaxy. 

Here we performed cosmological simulations to analyse the ef- 
fects of the star formation process inside stellar-gaseous disks 
embedded in a DM cosmological halo evolving in a fully con- 
sis tent scenario. We employ t he sub-grid star formation model 
bv lSpringel & Hernquistl d2003l) which is able to produce a self- 
regulated star formation in galactic disks. We point out that our 
model is not a general galaxy evolution model, since the grad- 
ual formation and growth of the stellar disk has not been taken 
into accou nt. However our approa ch, developed in two previ- 
ous papers dCurir et alj |2006. 2007), allows us to investigate the 
effects of different parameters like the disk-to-halo mass ratio 
and the gas fraction inside the disk, on the growth of the bar 
instability, its coupling with the star formation rate and its de- 
pendence on such parameters in a self-consistent cosmological 
framework. We compare the results of such a new set of simu- 
lations with th ose of our prev ious sets with the sta r formation 
switch ed off dCurir et alj ( 120061) . hereafter Paper 1, ICurir et ail 
d2007l) . hereafter Paper 2). In Paper 1 we presented simulations 
of purely stellar disks with the same disk-to-halo mass ratios, in 
the same cosmological scenario, and with the same initial condi- 
tions as in this new paper. In Paper 2 we investigated the growth 
and the evolution of the bar instability in disks with the same 
disk-to-halo mass ratios as in Paper 1, and different gas frac- 
tions, without star formation. Thus, we will compare the results 
of simulations here performed, with those corresponding to the 
same disk-to-halo mass ratio as in Paper 1 and to the same gas 
fraction as in Paper 2. The comparison between the parameters 
characterising the bar can be done only evaluating that of the old 
stars component, since the new stars are not present in Paper 1 
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and 2. In both such Papers it was shown that in DM dominated 
disks a bar feature is triggered and maintained by the Cosmology 
whereas, in the more massive disks, a gas fraction 0.2 is able to 
destroy the bar. The focus of the present work is to determine if 
the star formation changes such a result. 

The plan of the paper is the following: in Section 2, we sum- 
marise our recipe for the initial disk+halo system fully described 
in Paper 1, and present our star formation recipe. In Section 3 
there are our cosmological simulations, in Section 4 we point 
out our results. The parameters related to the bars formed in the 
new and the old stars and to the global, old+new, stellar popu- 
lations are given in such a section. Section 5 is devoted to our 
discussion and conclusions. It contains a Table including the pa- 
rameters characterising the old stellar bar at z= in the three 
Papers to allow a more easy comparison between the results of 
our Papers (Table 2). 

2. Method 

We embed a gaseous and stellar disk inside a cosmological halo 
selected in a suitable slice of Universe and follow its evolu- 
tion inside a cosmological framework: a ACDM model with 
Q m =0.3, Q.\ =0.7, <xg =0.9, h =0.7, where £2 m is the total matter 
of the Universe, £2a the cosmological constant, cr% the normal- 
isation of the power spectrum, and h the value of the Hubble 
constant in units of lOO/i -1 kms^Mpc -1 . 

A detailed description of our method to produce the cos- 
mological scenario where the disk is evolved has been given in 
Paper 1 . Here we present a short overview of our recipe. Paper 1 
showed that the numerical resolution does not impair their main 
result: in pure stellar disks, long living bars are a 'natural' out- 
come of the cosmological scenario. 

2.1. The DM halo 

To select the DM halo, we perform a low-resolution simulation 
of a concordance A CDM cosmological model, starting from 
redshift 20. 

With a standard 'friends of friends algorithm' we selected 
one suitable DM halo with a mass M-IO 11 /^ 1 M Q (at z=0). 
We resampl e it w ith the multi-mass technique described in 
iKlvpin et all (|2001|) . The particles of the DM halo, and those 
belonging to a sphere with a radius Ah~ x Mpc, are followed to 
their Lagrangian position and re-sampled to an equivalent res- 
olution of 1024 3 particles. The total number of DM particles 
in the high resolution corresponds to a DM mass resolution of 
1.21 10 6 /r' M Q . The high-resolution DM halo is followed to the 
redshift z=0. We run the DM simulation, to extract the halo 
properties in absence of any embedded stellar disk. The mass 
of our halo at z=0, 1.03 • 10 ll h~ l M , corresponds to a radius, 
Rvir = 94.7/z -1 Kpc, which entails 84720 halo particles. The se- 
lected halo is living in an under-dense environment. From its ac- 
cretion history (see Fig. 1 in Paper 1) we conclude that our halo 
undergoes no significant merger during the time it hosts our disk, 
nor immediately before. The halo density profile is well-fitted b y 
a Navarro, Frenk and White (NFW) form dNavarro et alj dl996): 
iNavarro et alj (1 19971) ) at z< 2. The concentration, C,,,>, equal to 
R v , r /iJ s , takes an high value, 18.1, c onfirming that this halo does 
"form" at quite high redshift ( e.g. IWechsler et alj d2002l) for a 
discussion about the link between concentration and assembly 
history of the halo). The dimensionless spin parameter of the 
halo is 0.04 at z=2, near to the average one for our cosmological 

model dMaller et al.ll2002tlBett et al]|2007l) " 

Cosmological DM haloes are not spherical dBett et alj2007l and 



references therein). From a rec ent analysis of th e halo shapes 
from the Millennium simulation. iBett et all d2007l) show that ha- 
los are triaxial and prolate. Their shape is driven by the hier- 
archical formation. The prolateness of our halo at z=2, where 
R v i r =30 Kpc, is 0.9, quite the same as at the disk radius (Paper 
1, Table 1). 

2.2. The baryonic disk 

The spatial distribution of particles follows the exponential sur- 
face density law: p = po exp -(r/r ) where ro is the disk scale 
length, ro = 4ft -1 Kpc, and po is the surface central den- 
sity. The disk is truncated at five scale lengths with a radius: 
Rrfijyt -20hT l Kpc. To obtain each disk particle's position ac- 
cording to the assumed densi ty distribution, we used the rejec- 
tion method dPress et alj|1986l) . We used 56000 star particles 
and 56000 gaseous particles to describe our disk. The (Plummer 
equivalent) softening length, the same for DM, gas, and star par- 
ticles, is 0.5/i~' Kpc in comoving coordinates. 

We embed the disk in the high resolution cosmological sim- 
ulation, at redshift 2, in a plane perpendicular to the angular mo- 
mentum vector of the halo and in gravitational equilibrium with 
the potential. Its centre of mass corresponds to the minimum po- 
tential well of the DM halo. The initial redshift corresponds to 
10.24 Gyr down to z = in our chosen cosmology. During the 
evolution, new star particles are formed from the gaseous parti- 
cles. We will refer to such new component as 'new stars compo- 
nent', whereas we will call 'old stars component' the non dissi- 
pative particles present in the disks at z=2. 

2.3. Star formation recipe 

We use the sub-grid star formation prescription by Springel & 
Hernquist (2003). In such prescription, when a gas particle over- 
comes a given density threshold its gas content is considered 
to reside in a multi-phase state, corresponding to the equilib- 
rium solution of an analytical model describing the physics of 
the multi-phase interstellar medium. Such a solution gives an ef- 
fective temperature for the gas particle, obtained as a weighted 
average of the (fixed) temperature of the cold phase and that of 
the hot phase. The temperature of the hot phase is set by su- 
pernova feedback and by the efficiency of cloud evaporation. 
The resulting inter stellar medium (ISM) has an effective equa- 
tion of state Pejjip) which is stiffer than isothermal and pre- 
vent Toomre instabilities even when a large amount of gas is 
present. When this prescription is used, the star formation rate 
agrees with the Schmidt law in disk galaxies, as obtained e.g. 
by Kennicut (1998). A detailed description of the numerical pre- 
scription and several tests can be found in Springel & Hernquist 
(2003). The effective temperature drives hydrodynamical inter- 
actions of the gas particles. Moreover, the model consistently 
gives a star formation rate, which is used to spawn a star par- 
ticle from the gas one on the basis of a stochastic prescription; 
the initial mass function adopted is the Salpeter's one. The star 
forming from the gas has the same position and velocity of the 
gaseous particle. 

3. Simulations 

We performed six cosmological simulations of a disk+halo 
system. We exploited the parallel Tree+SPH N-body code 
GADGET-2 dSpringelll2005l) (courtesy of V. Springel). The sim- 
ulations run on the CLX computers located at the CINECA com- 
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Table 1. Simulations: final values. o.s.=old stars , n.s.= new stars 



3 



N 


M disk 


gas fraction 


e(o.s.) 


e (n.s.) 


Q„ (o.s.) 


Qi (o.s.+n.s.) 


a„ ra _ v (o.s.) 


a,,,,,., (n.s.) 


bulge(o.s.) 


bulge(n.s.) 


bars in bars 


cl 


0.33 


0.1 


0.65 


0.72 


0.45 


0.5 
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c2 


0.33 


0.2 


0.55 


0.55 


0.51 


0.66 


11 


5.7 
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c3 


0.33 


0.4 


0.6 


0.55 


0.58 


0.69 


8.4 
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c4 


0.1 


0.1 


0.39 


0.01 


0.5 


0.41 
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c5 


0.1 


0.2 


0.45 


0.03 


0.46 


0.42 
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c6 


0.1 


0.6 


0.5 


0.02 


0.51 


0.47 
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puting center (BO, Italy ) and on OATo Beowulf-class cluster of 
32 Linux-based PC at the Osservatorio Astronomico di Torino. 
The main parameters and the final properties of our set of sim- 
ulations are listed in Table QJ We define as a geometrical mea- 
sure of the bar strength the value of the ellipticity, e = 1 - bja 
(Table [TJ; a strong bar corresponds to e > 0.4. A more dynami- 
cal measure of the bar strength at radius R has been defined by 

ICombes & Sanders] d 198 lb by using the parameter: Q, = < ^ (fi)> 
where Fj ax {R) = [d®(R, 6)/d8] max is the maximum amplitude 
of tangential force and < Fr(R) >- R(d<&o/dr) is the mean ax- 
isymmetric radial force derived from the m = component of 
the gravitational potential at the same radius, R. The maximum 
value of Q t (R) provides a measure of the bar strength Qt, for the 
whole galaxy. Stronger bars correspond to higher Qh values. In 
order to compare the results with the Papers 1 and 2 the values of 
Qi, are evaluated in the bars formed in the old stars component, 
since in our previous works we didn't include the star formation 
process. 

We give here in Table QJ also the value of the parameters re- 
lated to the bars formed in the new stars component and to the 
bar strenght of the global old+new stars populations. In TableQ] 
we present: the simulation number (I column), the mass of the 
disk in code units (i. e. 5.9 x 10 10 MO) (II col.), the fraction of 
gas, i.e. gas-to-disk mass ratio (III col.) Moreover we present, 
as final values (i. e. at z = 0) the maximum ellipticity for the 
old stars system(IV col.), and for the new stars (V col.) the bar 
strength according to Combe s~& Sandersl d 1 98 lb for the old stars 
(VI col.), and for the old+new stars system (VII col.), the major 
axis (physical Kpc) corresponding to the maximum ellipticity 
for the old stars (VIII col.)and for the new stars (IX col.), the 
morphology of the inner region of the disk for the old stars (X 
col.) and for the new ones (XI col.), peculiar features inside the 
disk (XII col.) 



4. Results 

In all our cosmological simulations, a stellar bar is still living at 
z=0 in the old star component. The new stars component at z=0 
is arranged in a bulge component which can or cannot present a 
barred shape depending on the initial mass of the disk and on 
the gaseous fraction. 

In Figs 1-8 we present the isodensity contours of the old stars, 
gas and new stars for simulations in Table QJ All such figures 
have been built up with the same box-size, number of levels, 
and density contrast (see caption of Fig.QJi, as in Papers 1 and 2. 



4.1. Morphologies 

We will discuss first the bar features as far as the old star com- 
ponent is concerned. FigQJand Figf2]show the evolution of this 
component from z=2 to z=0 for simulation cl and c4. The bar 
feature is well defined since z=l .75 in both the cases. In the more 



















Fig. 1. Face-on isodensity contours of simulation cl from z=2 
to z=0. Spatial resolution is always 0.5h~ l physical Kpc and 
the box size is 40 times the spatial resolution. Contours are com- 
puted at 1 1 fixed levels ranging from 2 x 10~ 4 to 0.015 in term 
of density fraction of stars (gas) within the spatial resolution to 
the total star (gas) density in the map. 



massive case, the evolving bar becomes stronger and dominates 
the disk, whereas in the less massive case, the bar remains a 
central structure inside the disk, generating a bar-in-bar feature. 
The strongest bar arises from simulation cl (Fig. [3). In this case, 
the ellipticity, evaluated through the isodensity plots, is higher 
(0.65) than the one measured in the purely stellar disk of the 
same mass (0.52, Paper 1), but slightly smaller than that mea- 
sured in the case of a disk with the same mass and the same 
gas fraction (0.68, Paper 2). The enhancement of the stellar bar 
strength in comparison to the gas free case is due to the coupling 
of the gaseous component to the stellar one: the gaseous bar is 
indeed superimposed to the stellar one, as pointed out in Paper 
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Fig. 2. Face-on isodensity contours of simulation c4 from z=2 
to z=0. Isodensity levels and spatial resolution are the same as 
in Fig.Q] 



2. However, this effect is slightly reduced by the star formation 
activity. 

In the simulation c2, the stars produce a barred feature, last- 
ing until z — (Fig. [4]). This is an important difference with the 
case having the same gas fraction in Paper 2, where the bar was 
decreasing its strength and disappearing at z ~ 0.15. The same 
behaviour arises also in the c3 case (Fig. [5]), which corresponds 
to a larger initial gas fraction: the bar is living until the final 
redshift, whereas, without the star formation, the bar was disap- 
pearing at z ~ 0.6. 

The gas forms a bar since the beginning of the disk evolution and 
the new stars originating from the gaseous component assume an 
elongated structure which fuels the old stars bar feature: the new 
stars are dynamically cold and they strongly feel the resonances 
in the equatorial plane, where they are born. In the simulations 
with higher gas fractions, they are swept out in the central re- 
gions by vertical reson ances which results i n the formation of 
a central young bul ge (Fried li & Benz|[l995l) . According to the 
claim of iBerentzen et al.1 d 1998b . the growth of a gas mass con- 
centration in the disk dissolves the regular orbits in the stellar 
bar. In Paper 2, we deduced a threshold value for the gas con- 
centration inside a radius of 2 Kpc, 9%, able to destroy the bar 
in the more massive disks. If the star formation is included, it 
transforms the dissipative component into a non dissipative one 
and therefore the central concentration of the gas is lower. Fig. 



Fig. 3. Face-on, edge-on and side-on isodensity contours (from 
left to right) of simulation cl at z=0 (see text); top panel shows 
the gaseous component, middle panel the old star component, 
bottom panel the new formed stars. Isodensity levels and spatial 
resolution for gas and stars are as in Fig.Q] 
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Fig. 4. Same as in Fig.[3]but for simulation c2 
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Fig. 5. Same as in Fig.[3]but for simulation c3 
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Fig. 6. Same as in Fig.[3]but for simulation c4 



Fig. 8. As in Fig.[3]for simulation c6. 




Fig. 7. As in Fig.[3]for simulation c5. 



|9]and Fig.fTOlshow the evolution of the gas mass in the inner re- 
gions, i.e., of the gas concentration (upper panel), and that of the 
mass of gas plus new stars (lower panel) in the same regions, in 
simulations c2 and c3 respectively. In both these cases, the gas 
mass inside 2 kpc never overcomes the threshold value pointed 
out in Paper 2. However, the gas+new star concentration in sim- 
ulation c3 overcomes this threshold value at z=0.6, without any 
consequence on the bar observed (Fig. ITTb . This is because the 
effect of bar dissolution due to the central mass concentration 
cannot be achieved including the mass of new stars which are 
crossing the central region through elongated orbits and are not 
really concentrated as the gas. In the simulation c3 the gas shows 
more developed arms at z=l than in the case cl (Fig.ll21i. But 
this feature, which is considered as a possible mechanism for the 
bar destruction has n o impact on the bar in these new simulations 
dBournaud et al.ll2005l) . 

A different discussion is devoted to the DM-dominated disks, 
which, as we already stressed in Paper 1 and 2, show bars which 
are not a classical product of the self-gravity or of angular mo- 
mentum exchanges, but they are features that strongly depend on 
the dynamical state and evolution of the cosmological halo, and 
we did not observe their distruction even by an higher gas con- 
centration. Old stars in simulation c4 show a stellar bar stronger 




redshift 



Fig. 9. Gas-to- disk mass ratio at different disk radii : 2 Kpc 
(magenta dashed-dotted line), 3 Kpc (blue dashed line), 4 Kpc 
(red dotted line), 5 Kpc (black full line) for simulation c2. In the 
upper panel we present the time behaviour of the gas as fraction 
of the disk mass, in the lower panel we see the same for the sum 
of the gas and the new formed stars. 



than that in the pure stellar disk of the same mass (ellipticity 
0.39 instead of 0.3, see Tab. [2]) but weaker than that arising 
from the simulation with the same disk mass and gas fraction but 
without star formation (Paper 2). Bar-in-bar features are present 
along the whole evolution (simulation c4) until z=0 (Fig. [6j. 
Simulations c5 and c6, which correspond to increasing gas frac- 
tions, show similar features during the system evolution. In this 
set of DM dominated disk simulations, the final bar is increasing 
its strength with the increasing gas fraction (Fig.Q~3] Tab.[TJi. This 
trend is different from the one we observed in Paper 2, where the 
increasing gas fraction produces an increase to a maximum and 
then a decrease of the bar strength. 

In DM dominated disks the bars in the old star component are 
the imprint of the halo cosmological properties, as discussed in 
the previous Papers. Fig. [14] shows the isodensity contours of 
the halo in the plane of disk in the same frame, with the same 
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2 1.5 1 0.5 

redshift 



Fig. 10. Same as in Fig.[9]but for simulation c3. 




2 1.5 1 0.5 

redshift 



Fig. 11. Evolution in redshift of the ellipticity for simulations c2 
(blue full line) and c3 (red dashed line). 

isodensity levels, and box sizes as those of the old stars in Fig. 
Q] and Figj2] The semimajor axis of the triaxial DM halo, for 
Z < 1 .8 is both in phase with the central bar of the disk and they 
have the same pattern speed. 

In the more massive disks (Fig. [T5| ), the central DM region 
becomes more round, due to impact of the disk itself, and the 
coupling with the stellar bar disappears. Therefore, in the case 
of the lighter disks, the cosmological halo properties drive the 
formation of the bar, whereas in the more massive disks, the bar 
is a classical product of the disk instability. In the DM dominated 
disks the new stellar component produces a small spheroidal 
bulge which has a disky shape in the simulation c6. The gas com- 
ponent is confined in the central regions from the beginning of 
the evolution, and therefore the new stellar component arranges 
itself in a small spheroidal bulge. Therefore, in DM dominated 
disk, the new stars at z=0 have no bar features, at difference with 
more massive disks. We recall that, in DM dominated disks the 
bars are not due to classical resonances. So, the new stars do not 
feel these resonances whereas this is what happens in the cases 




Fig. 12. Gas morphologies at z=l for simulations cl (left panel) 
and c3 (right panel). 



0.8 




0.2 r 



nn l i i i 

0-0 0.' 0.2 0-3 0.4 



Fig. 13. Behaviour of the bar strength and ellipticity at z= for 
our set of cosmological simulations with increasing gas fraction. 
Qb (dotted line) and ellipticity (full line) of our more massive 
disks (i.e. disk-to-halo mass ratio 0.33); Qb (dot-dashed line) 
and ellipticity (dashed line) of our less massive, DM-dominated, 
disks (i.e. disk-to-halo mass ratio 0.1). 



of the more massive disks. The decoupling between new and old 
stars is depicted also by In Fig. [I7]the distribution of radial ve- 
locity dispersions in a massive disk at z=l is compared with the 
one in a DM dominated disk at the same redshift; in such a case 
the decoupling between new and old stars is well depicted. The 
radial velocity dispersion of this new star component is lower 
than that of the old stars. At the end of the evolution we there- 
fore see a small bulge of new stars hiding a small nuclear bar of 
old stars. Such a nuclear bar remains, as observed in Papers 1 
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Fig. 14. Isodensity contours of the DM in the xy plane, from 
redshift z = 2 to z = for the simulation c4. Density levels and 
spatial resolutions are the same as in Fig. Q] The syncronicity 
of the semimajor axis of the halo with the one of the inner old 
stellar bar can be noticed by direct comparison of each panel 
(after z=1.8) with the corresponding one in Fig. [2] 



and 2 and stressed by Fig. [14] as a genuine product of the cos- 
mology, since the light disks are classically below the instability 
threshold. 

4.2. Bar pattern speed 

The bar pattern speed, D. p is the angular speed of the bar-like 
density wave as viewed from an inertial frame. Here, it is evalu- 
ated following the position angle of the major bar axis during 
the disk evolution, where t he major axis is selected with th e 
same criterium El quoted in Miche l-Dansac & Wozniakl f2006). 
namely as the radius where the ellipticity profile reaches a max- 
imum as in Paper 2. The values of Q. p as function of the redshift 
for the more massive disks and the lighter disks are shown in 
Fig. [18] and Fig. [Irrespectively. 

For the more massive disks, we observe a decrease of Q p in 
all the simulations. This is a different behaviour from the one 
observed in Paper 2: there we noticed an increase of the pat- 
tern speed coupled with the increasing gas fraction, since the 
gas presence both shortens the bar and decreases its ellipticity. 



Fig. 15. Isodensity contours of the DM, as in the previous figure, 
for the simulation cl. 

This mechanism is no longer working when the star forma- 
tion is activated, since the ongoing star formation reduces the 
gas concentration in the central regions and the new stars ar- 
range into an elongated structure which enforces also the length 
of the original bar (Fig. [20] see also the previous section). 

For the DM dominated disks the evaluation of Q, p in this pa- 
per as in the previous ones, is done by referring to the angular 
velocity of the long-lived central bar. Its behaviour is the same as 
observed in Paper 2, namely it shows a fast decrease before red- 
shift z — 1 , and then a quasi stationary behaviour until redshift 
0. 

5. Discussion and conclusions 

We presented six cosmological simulations with the same disk- 
to-halo mass ratios as in Paper 1 and Paper 2. In order to study 
the impact of the forming stars from the gaseous component, 
here we included and varied its percentages inside disks of dif- 
ferent disk-to-halo mass ratios, as in Paper 2 where the star for- 
mation was switched off. In Table [2] we show the final data for 
ellipticities and semimajor axes resuming the results of our three 
papers. 

The old star component show a long lasting bar, 10 Gyr old, 
in all the simulations of this work. Moreover, in all such simula- 
tions, except for simulation c6, the bar is stronger than that de- 
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Table 2. Cosmological disks simulations: global results 



Disk mass 


gas fraction 


no star formation 
ellipticity 


@>max 


star formation 
ellipticity 




0.33 


0. 


0.52 


5 






0.33 


0.1 


0.68 


8.4 


0.65 


8 


0.33 


0.2 


0.1 


no bar 


0.55 


11 


0.33 


0.4 


0.07 


no bar 


0.6 


8.4 


0.1 


0. 


0.3 


6.5 






0.1 


0.1 


0.58 


5.8 


0.39 


3 


0.1 


0.2 


0.6 


5.4 


0.45 


3 


0.1 


0.6 


0.42 


5.8 


0.5 


3 





Fig. 16. Same as in Fig. Q~2]but for simulations c4 and c6. 



veloped in the pure stellar case with the same disk-to-halo mass 
ratio but weaker than that formed in the case of the same gas 
fraction without star formation. 

We find that the star formation, reducing the central gaseous 
mass concentration, allows the bar to survive until the end of the 
evolution also in the more massive disks, at variance with the re- 
sults in Paper 2 (where the gas was not allowed to form stars): in 
such a case a gas fraction 0.2 was able to destroy the bar. Even if 
some details of the disk morphologies obviously depend on the 
star formation prescription, the bars arising here in the c2 and c3 
simulations are due to the reduced central gas concentration and 
thus to the weakening of its effect on the bar itself. We verified 
this point by rerunning simulation c2 changing the star forma- 
tion prescription. Instead of the GADGET-2 effective model, we 
used a simple algorithm in which, when a gas particle reaches 
a given density threshold with a temperature lower than a min- 
imu m temperature thr eshold, it is converted into a star particle, 
as in lKatz et al.l d!996l) . Fig.l2Tlshows the morphology of the old 





Fig. 17. Radial velocity dispersion at z=l of old stars (red open 
circles) and new stars (green crosses) for simulations c2 (left 
panel) and c5 (right panel). 



stellar component for such a different star formation recipe: the 
survival of the bar is due to the star formation, independently on 
its details. 

In all the simulations of more massive disks, the new stellar com- 
ponent shows a barred shape coupled with that of old stars but 
with semimajor axis and ellipticity values smaller than those. In 
all the simulations of DM dominated disks, a bar feature is main- 
tained in the old stellar component at the disk center, and the new 
stars form a spheroidal bulge that hides the bar of the old stars. 
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.5 1 0.5 
redshift 



15 




1.5 



1 

redshift 



0.5 



Fig. 18. Evolution of the bar pattern speed with the redshift for 
simulation cl (black dashed line, c2 (red full line) and c3 (ma- 
genta dotted line). 



Fig. 20. Behaviour of the semimajor axis for the simulations c2 
(blue full line) and c3 (red dashed line). 




.5 1 0.5 

redshift 



Fig. 19. Evolution of the bar pattern speed with the redshift for 
simulation c4 (black dashed line), c5 (red full line) and c6 (ma- 
genta dotted line). 



This effect could be of interest as far as the observations about 
the Milky Way are concerned where, in addition to the main bar, 
a small nuclear bar c onsisting of o ld stars seems to be included 
in the bulge (see e. g. lAlardl ((20011) ). 

The final strenght of these old star bars increases by increas- 
ing gas fraction and their pattern speed is quickly decreasing 
before z= 1 . 

The classical results obtained outside the cosmological sce- 
nario are no longer applicable. This conclusion remarks the re- 
sults of Paper 1, where it was shown that in the DM dominated 
disks the bar feature is triggered and maintained by the cosmo- 
logical properties of the halo (namely its triaxiality and its dy- 
namical state). 

From the models presented here we suggest that a very low 
pattern speed (few Km s~ l Kpc~ l ) could be a signature of a dom- 
inating halo and not a classical product of the disk instability. 




Fig. 21. Isodensity contours of the old stars at z=0, for the sim- 
ulation c2, rerunned with the Katz'96 star formation recipe 



The whole set of cosmological simulations we presented 
here is suggesting that the star formation works in favour of 
maintaining the bar feature in self gravitating disks, against the 
effect of increasing gas fractions, whereas in DM dominated 
disks, it slightly reduces the bar strength. Moreover, inside such 
disks, the new stars are arranged in non barred systems which 
are decoupled from the bar of the old star component. 
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